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Abstract.  The  accuracy  of  the  rigid-Earth  solution  SMART97  is  2  //as 
over  the  time  interval  (1968,  2023).  To  obtain  a  nonrigid-Earth  solution, 
we  use  the  transfer  function  of  Mathews  (1999).  The  perturbations  of  the 
third  component  of  the  angular  velocity  vector  are  taken  into  account. 


1.  The  rigid-Earth  solution  SMART97 

The  differential  equations  of  the  rigid-Earth  rotation  are 
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The  torque  (L,  M,  N )  has  to  be  computed  by  using  simultaneously  the  zonal 
and  the  tesseral  harmonics  because  the  first  and  the  second  derivatives  of  the 
diurnal  and  the  semidiurnal  terms  are  very  important  with  respect  to  the  long- 
period  terms.  Table  1  gives  the  amplitudes  of  the  semidiurnal  term  (coming 
from  0*2,2  an(i  ^2,2)7  of  the  18.6-year  term  and  of  the  13.66-day  term. 

Integrated  in  this  way,  the  rigid  Earth  solution  SMART97  (Bretagnon  et 
al.,  1998)  can  reach  a  high  accuracy.  It  has  been  compared  with  a  numerical 
integration  using  DE403  (Standish  et  al.,  1995)  for  the  positions  of  the  Moon, 
the  Sun  and  the  planets.  The  accuracy  is  2  //as  over  1968-2023.  Figure  1  gives 
the  differences  for  ip,  lo  and  for  the  Earth  rotation  angle  Lp. 
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Figure  1.  Theory  SMART97  —  numerical  integration  using  DE403 
over  1968-2023. 


2.  The  influence  of  the  truncation  level 

Table  2  gives  the  number  of  the  periodic  terms  for  a  =  O.Ol^as,  a  =  0.1/uas,  and 
a  =  l^as  and  the  number  of  Poisson  terms  greater  than  a  over  [J2000. 0  —  100  yrs, 
J2000. 0+100  yrs],  for  the  rigid-Earth  solution  SMART97.  The  table  also  gives 
the  accuracy  of  the  solution  for  different  levels  of  truncation.  Data  are  given 
for  the  variable  p  (At/’  in  the  IERS  Conventions  1996)  of  nutation  in  longitude 
reckoned  from  the  equinox  of  date.  The  number  of  terms  of  the  Herring  solution 
for  a  nonrigid  Earth  (McCarthy,  1996)  is  given  by  comparison. 


Table  1.  Amplitude  of  the  18.6-year,  18.66-day,  12-liour  terms  and 
of  their  first  and  second  derivatives. 

period  (/’(in")  (/’(in  "/yr)  (/’(in"/yr2) 

18.6  years  17.280776  5^38  OfT 

13.66  days  0.221507  37.212  6251.37 

12  hours  0.000036  0.132  762.56 


The  Poisson  terms  In  SMART97,  the  Pxt  Poisson  terms  (Poisson  terms  of 
degree  1)  represent  13.6%  of  the  periodic  terms  with  a  truncation  of  l^as,  16.8%) 
with  a  truncation  of  0.1/uas  and  20.8%)  with  a  truncation  of  0.01/uas.  That  must 
be  compared  to  the  Herring  solution  in  which  the  Pxt  Poisson  terms  represent 
29.1%)  of  the  periodic  terms  with  a  truncation  of  lpas. 
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The  truncation  In  the  construction  of  the  SMART  solution  for  a  nonrigid  Earth 
using  the  transfer  function  of  Mathews,  the  amplitudes  of  some  terms  in  reso¬ 
nance  can  be  increased  considerably.  So,  to  obtain  all  the  terms  greater  than  0.1 
gas  in  the  nonrigid-Earth  solution  it  is  necessary  to  keep  a  level  of  truncation 
of  0.01  gas  for  the  rigid-Earth  solution.  We  give  the  number  of  terms  of  the 
SMART  solution  for  a  nonrigid  Earth  in  Table  3. 


Table  2.  Number  of  terms  of  the  nutation  in  longitude  p  in  SMART97 
for  different  truncation  levels  a. 


Solution  and  a 

P 

Pxt 

Pxt2 

hd 

X 

c-4- 

co 

Pxt4 

Total 

Accuracy 

SMART  (0.01/ias) 

3910 

815 

183 

14 

2 

4924 

2.2  gas 

SMART  (0.10/ias) 

1586 

266 

61 

4 

0 

1917 

8  gas 

SMART  (1.00/ias) 

642 

87 

17 

1 

0 

747 

40  gas 

Herring  (1.00/ias) 

375 

109 

0 

0 

0 

484 

Table  3.  Number  of  terms  of  the  nutation  in  longitude  of  the  nonrigid 
Earth  solution  at  the  0.1  gas  truncation  level. 

Periodic  Pxt  Pxt2  Pxt3  Pxt4  Total 

1581  264  60  4  0  1909 


3.  Nonrigid-Earth  solution 


3.1.  The  transfer  function  of  Mathews  (1999) 

To  compute  the  nonrigid-Earth  solution,  we  use  the  transfer  functions  in  a  strict 
process.  For  instance,  the  transfer  function  of  Mathews  (1999)  is 


T(<J-  e\eR) 
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with  eR  =  0.003284507,  (1  +  a)  is  the  nutation  frequency  in  space.  We  use  the 
complex  frequencies 

(T i  =  (  0.0025940,  -0.0001438546),  a2  =  (-1.00231861,  0.00002578), 
as  =  (-0.998957,  0.000687),  a4  =  (  0.000413499,  0.000000280), 


and  the  complex  coefficients 

Ao  =  (  1.0000099,  -0.37652854  x  10“®), 

IVi  =  (-0.79952969,  0.043796154), 

N2  =  (  0.048964919,  0.16332679 X  10“2), 

N3  =  (  0.29445472  x  10“3,  -0.82328898  X  10“4), 
N4  =  (-0.15139223  X  10“4,  -0.11248592  x  10“5). 
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3.2.  Application  to  the  Earth’s  angular  velocity  vector 

We  apply  the  new  transfer  function  not  to  the  quantities  sin  eoA ip  and  Ae  but 
to  the  Earth’s  angular  velocity  vector  (p,  q,  r)  of  the  rigid  case  expressed  as  a 


function  of  the  three  Euler 

angles 

PR  = 

ipR  sin  uR  sin  ipR  +  uR  cos  ipR, 

QR  = 

ipR  sin  uR  cos  ipR  -  uR  sin  ipR, 

(1) 

tr  = 

IpR  COS  UR  +  ip R  . 

For  each  prograde  or  retrograde  argument,  the  angular  velocity  vector  of  the 
nonrigid  Earth  is  obtained  by  multiplication  by  the  transfer  function 

(pnr  +  i  QNR )  =  (pr  +  i  Qr)  x  T{cr;  e\eR). 

From  the  quantities  ( Pnr ,  qNRi  rNR ),  we  can  compute  the  derivatives  of  ip,  u 
and  ip  in  the  nonrigid  case  by  the  inverse  of  (1) 

i’NR  =  ( PNR  sin  ipNR  +  qNR  cos  (pNR)/  sin  uNR , 

d^NR  =  (pnr  cos  ip]\fR  Qnr  sin  pnr),  (2) 

Pnr  =  rNR  -  ipNRcosuNR. 

These  equations  are  strict  but  the  right-hand  members  depend  on  u  and  <p  and 
the  derivative  of  ip,  and  we  have  to  proceed  by  iterations  to  solve  this  system. 
The  process  converges  without  difficulties.  The  precision  of  the  computation  of 
ipNR  is  better  than  3.5  X  10_9///year  which  yields  an  accuracy  of  0.01  pas  for  the 
18.6-year  term  and  0.50  pas  for  the  883  year  term.  The  classical  method  is  to 
assume  tnr  =  rR.  The  results  are  illustrated  by  the  following  table.  We  denote 
SM97M99  as  the  solution  obtained  by  SMART97  +  Mathews  (1999).  We  can 
see  the  results  from  Mathews’s  function  are  very  close  to  the  series  of  Herring 
(McCarthy,  1996). 


Solution 

Argument 

p  (sin) 

p  (cos) 

e  (sin) 

e  (cos) 

SM97M99 

A3  +  D  —  F 

17  206  664 

-3  357 

-1488 

-9  205  156 

Herring 

17  206  394 

-3  702 

-1523 

-9  205  474 

SM97M99 

2  A3 

-1318  625 

-662 

-486 

573 040 

Herring 

-1318  526 

-670 

-471 

573 046 

SM97M99 

2A3  +  2  D 

-227  663 

309 

150 

97  854 

Herring 

-227  720 

269 

136 

97  864 

SM97M99 

a3 

-36  674 

-123  068 

16  616 

684 

Herring 

-36  777 

-123010 

16  590 

698 

The  full  method  is  to  take  into  account  the  tidal  variations  in  the  Earth’s 
rotation.  We  introduced  the  series  (u  —  us)  and  (u  —  (McCarthy,  1996)  in 
order  to  obtain  r^R  from  rR, 


rNR  =  rR  +  (oj  -  Us)  +  (u  -  UD). 
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Figure  2.  Influences  of  (uj  —  ujg )  and  (uj  —  ago  )  on  p  and  e  over  1968- 
2023.  Units  are  pas. 


The  modifications  are  very  important  not  only  for  the  third  Euler  angle  <p  but 
also  for  the  first  two.  For  instance,  we  obtain 

A p  =  42  242  /  +  69  sin(2A3  +  2D  -  2F)  -  4  sin(  A3  +  D  -  F)  +  ■  ■  • , 

Ae  =  — 15  /  +  712  cos(  \3  +  D  —  F)  —  19  cos(2A3  +  2D  —  2F)  +  •  •  • , 

where  the  amplitudes  are  in  pas  and  the  time  in  thousands  of  Julian  years  from 
J2000.  Therefore,  for  s,  the  18.6-year  term  becomes 

e  =  —  1  488  sin(  A3  +  _D  —  F)  —  9  204  444  cos(  A3  +  D  —  F ) . 

By  comparison  with  the  previous  table  we  see  the  necessity  to  determine  again 
the  transfer  function.  The  modifications  of  p  and  e  are  plotted  in  Figure  2. 


4.  Conclusion 

We  have  now  to  determine  a  new  nonrigid-Earth  rotation  solution  by  using  the 
new  model  of  Mathews  et  al.  (2000)  taking  into  account  the  tidal  variations  in 
tli  Earth’s  rotation. 
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